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Tij- . Abstract 

0^ 

■^ I The production of jet pairs with small transverse momentum and large rel- 

o ■ ative rapidity in high energy hadron-hadron collisions is studied. The rise of 

JL , the parton-level cross section with increasing rapidity gap is a fundamental 

q_) I prediction of the BFKL 'perturbative pomeron' equation of Quantum Chro- 

1-^ ' modynamics. However, at fixed collider energy it is difficult to disentangle this 

effect from variations in the cross section due to the parton distributions. It 
k><( \ is proposed to study instead the distribution in the azimuthal angle difference 

;h ' of the jets as a function of the rapidity gap. The flattening of this distribution 

with increasing dijet rapidity gap is shown to be a characteristic feature of the 
BFKL behaviour. Predictions for the Fermilab pp collider are presented. 



1 Introduction 

There is currently much interest in the QCD 'perturbative pomeron'. This is the 
phenomenon, obtained by resumming a certain type of soft gluon emission to all or- 
ders in the leading logarithm approximation using the Balitsky-Kuraev-Fadin-Lipatov 
(BFKL) equation ^, which is supposed to produce a sharp rise at small x in deep 
inelastic structure functions 0. Recent measurements at HERA may well show 
the first evidence for this behaviour, although it is premature to draw any definitive 
conclusions. 

One of the difficulties with extracting information on the perturbative pomeron 
from structure function measurements alone is that both perturbative and non- 
perturbative effects are very likely intertwined in a non-trivial way, and therefore 
obtaining information on the former requires some sort of model-dependent subtrac- 
tion of the latter. For this reason, attempts have been made to find other quantities 
which probe more directly the perturbative behaviour. In the context of deep inelas- 
tic scattering, one can look for an associated jet of longitudinal momentum fraction 
x' ^ Xbj and comparable transverse momentum to the virtual photon momentum Q 
^. The BFKL behaviour is then refiected in the growth of the cross section with 
log(x'/xBj). Alternatively, Mueller and Navelet have shown that in high-energy 
hadron-hadron collisions one can utilize a two-jet inclusive cross section where the 
jets have small transverse momentum and a large relative rapidity. Ay. The rise in 
the cross section with increasing Ay is then controlled by the perturbative pomeron. 

It is this latter process that we investigate here. We first define the cross section 
introduced in Ref. [^, and show why it is difficult to measure under present experi- 
mental conditions. We then consider an alternative but closely-related quantity, the 
azimuthal angle correlation between jets at large relative rapidity, which we argue is 
more promising from an experimental point of view. Predictions are presented for the 
Fermilab pp collider at ^i = 1.8 TeV. Our results also apply to the corresponding 
two-jet cross section in deep inelastic scattering as measured at HERA, although we 
shall not pursue this issue here. 

Consider the inclusive two-jet cross section in pp (or pp) collisions at energy ^/s, 
where each jet has a minimum transverse momentum M -C i/i, and the jets are 
produced with equal and opposite large rapidity ±|A.0 Adapting the results of 
Ref. IQ, the cross section for this can be written in the form 



da 



dyidy2 



yi=-y2=2^ 



xiG(xi, M^) X2G{x2, M^) a{asiM^), M\ A) 



^The restriction to equal and opposite rapidities is not crucial and is made here only to simplify 
the discussion. 



where 



4 
G(a;,/i^) = gix,fx^) + - ^(g(x,/i^) + g(x,/i^)) 

9 q 

a;i = X2 = ^ cosh(iA) ~ —=e^', 
's Js 



a{as,M\A) 



a.CA\' n^ 



TT / 2M2 



1 + ^ a„ (a,A)" + . . . 



n>l 



(2) 



Note the use of the 'effective subprocess approximation', appropriate here because of 
the dominance of small momentum transfer in the subprocess t channel. The ... in 
(0) refers to corrections outside the leading logarithm approximation implicit in ([^), 
i.e. terms of order q;"A"^^, a"A"^^, . . . and power-correction terms suppressed by 
powers of e~^. 

According to the BFKL 'perturbative pomeron' analysis, the subprocess cross 
section M^a is expected to rise at large A. The asymptotic behaviour, ignoring 
possible 'parton saturation' effects 0, is predicted to be 

M^a ~ e^^, as A ^ oo , (3) 

with A a number of order 0.5. This is the analogue of the expected x"^ growth of the 
F2 structure function at small x 0. 

Now with ^/s = 1.8 TeV and M ~ 10 GeV, the maximum value of A is of 
order 10, which might at first sight appear sufficiently large to test the predicted 
BFKL behaviour.0 The problem, however, is the additional A dependence in the 
cross section induced by the x dependence of the parton distributions in (|l|). At large 
A, the behaviour of the cross section will be completely dominated by the x -^ 1 
suppression of the parton distributions. This can only be avoided by increasing the 
energy i/i of the collider as A is increased, in such a way that Xi and X2 remain fixed, 
which is not an easy proposition in practice. 

There are at least two ways around this difficulty. First, one could argue that the 
parton distributions at medium-to-large x are now sufficiently well known that they 
can be factored out of the measured cross section with sufficient accuracy. While this 
might be true for the quark component of G{x, fi"^), the gluon component is much less 
well known at large x. Even if it was, there is still a problem with scale dependence 
— the cross section in (|I|) is not known beyond leading logarithm accuracy, and so 
the choice of scale in the parton distributions is somewhat arbitrary. 

A second possibility is to use the distribution in the azimuthal angle difference 
between the two jets as a signature for BFKL behaviour. When A is small, the cross 

^In fact, the predicted asymptotic behaviour appears to set in at relatively modest values of A, 
see Section 2 and Ref. ||^. 



section is dominated by the lowest order 2^2 scattering subprocesses and the jets 
are back-to-back in the transverse plane. As A increases, more and more soft gluons 
with transverse momentum ~ M are emitted in the rapidity interval between the two 
fast jets, and the azimuthal correlation is gradually lost until, asymptotically, there 
is no correlation at all. The change in the overall weighting, due to variations in the 
parton distributions as A increases at fixed ^/s^ has no effect on the shape of the 
azimuthal distribution. 

The study of the azimuthal correlation between the two jets at large rapidity is 
the subject of the present analysis. We first of all derive the basic formulae, which 
is an extension of the treatment in |^, and then make numerical predictions for 
experimentally measurable quantities. We shall demonstrate that the weakening of 
the correlation should already be observable at the Fermilab pp collider. 

2 BFKL formalism for dijet production 

We start by considering the subprocess cross section for inclusive two-jet (i.e. par- 
ton) production. We are interested in jets produced with equal and opposite large 
longitudinal energy, k^ ^ E where E is the parton beam energy in the parton-parton 
centre-of-mass, and small transverse momentum kx > M where M ^ E' is a fixed 
cut-off. The rapidity gap between the jets is then Ay = A ~ 21og(ii^/M) ^ 1. In 
practice, M will be a number of order 10 GeV and the maximum value of E is set 
by the kinematic limit ■\/i/2. 

At the Fermilab collider, a large fraction of such jet pairs are produced in gluon- 
gluon collisions, and so in what follows we focus on the subprocess gg ^ gg + X. 
The quark initiated processes are reinstated afterwards using the effective subprocess 
approximation. If we assume, to begin with, that we can ignore the running of the 
strong coupling and take as = as{M'^) as fixed, then the subprocess cross section of 
Eq. (|lD can be written as @ 

aKM^A)=f^' -^ rd0F(0,A). (4) 



n J 2AP J-n 

Here we have introduced the variable = vr — (pjj where (pjj is the azimuthal angle 
difference between the two jets. Thus = corresponds to back-to-back jets in the 
transverse plane. In what follows we will be particularly interested in the differential 
distribution da/d(j), which is proportional to F. 

The function F in Eq. (^) is simply the quantity /(/cri, ^T25 ^) defined in Ref. 0, 
integrated over M^ < k'^^ < oo at fixed 0. The latter function satisfies the BFKL 
equation (see Appendix). The solution for F is 

F{<P,A) = - Y: e^"'^C7„(t), t = ^^A, (5) 



where 

1 /•+°° dz 



X„W = Re{i,(l)-iii,(l + \n\) + tz)], (6) 

and iIj{x) is the logarithmic derivative of the gamma function. Note that this re- 
suh corresponds to a perturbative expansion in powers of a^A ~ aslog^E^/M"^). 
It is these leading logarithms which have been resummed by the BFKL equation. 
Substituting back in the original expression for the cross section gives 



da 



dyidy2 



c^ xMxi,M')x2G{x2,M') (^) _^F(0,A), (7) 
yi=-y2=2^ 



with xi = X2 = 2Mcosh(|A)/y^. Before performing a full numerical calculation 
of this cross section, we discuss several important analytic results which follow from 
Eqs. (|,D. 

2.1 Comparison with exact lowest-order calculation 

Application of the BFKL formalism only makes sense in a kinematic region where 
the exact leading order (2 -^ 2) cross section is well-approximated by the first term 
in the perturbation series on the right-hand side of Eq. (0). The latter is obtained 
by setting t = in Eq. (pf), which gives 



Cn{0) 



1 /■+°° dz 



27r J-oo z^ + \ 



=^ F{<P,A)=6{<P). (8) 

When this is substituted in Eq. (|^), we obtain 

M' ^Ilo = Ivra^Cl . (9) 

Note that in the A ^ oo limit the subprocess cross section scales as 1/M^. In 
this section we rederive this result starting from the exact (2 -^ 2) lowest order 
cross section, which will enable us to assess how rapidly the asymptotic behaviour is 
approached. 

The two-jet inclusive cross section in leading order is given by 



da 1_ 

a,b,c,d=g,g 



dyidy2dpj. lovr^s 



where J2 denotes the appropriate sums and averages over colours and spins. To begin 



with, we restrict our attention to the gg -^ gg subprocess. Setting yi = —y2 
gives 

xg{x,iJ,^) xg{x,fi^) 



da 



dyidyidp 



yi=-y2=^^ 



25Qttp^ cosh^(iA) 



T.\M{gg^gg)\' 



where 



X 



2pj 



cosh(iA) 



and the subprocess matrix element is evaluated at 
1 



- = --(l-tanh(iA)) , - 

Integrating over the jet transverse momentum pt > M then gives 



--(l + tanh(iA 



da 



dyidy2 



yi=-y2- 



M2 Pj, 



2567r cosh^(iA) 



2m2 



iA 



'11^ 



(12) 



(13) 



(14) 



with p|.(max) = ;|scosh~ (|A). By extracting a factor of M~^, and ignoring loga- 
rithmic scaling violations in the parton distributions, the integral on the right-hand 
side becomes a function of the dimensionless quantity X = 2Mcosh(|A)/y^. 



da 



E\Migg^gg)\' 



dyidy2 

Next, consider the limit 



X-' du' 



, , 1. 25677 cosh^(iA) M2 A 



M" 



xg{x,iJ^) 



x=Xu 



K A < loe 



M . 



The first term on the right-hand side of Eq. (|15]) tends to a finite value of 

E\M{gg—'^' 



2567r cosh^(iA) 






(15) 



(16) 



(17) 



The integral over the parton distributions is dominated by the contribution from the 
lower limit, i.e. 



^-' du^ 



xg[x, fi 



x=Xu 



n u 

Combining the results from Eqs. (1^,18) gives 

da 



dl/id|/2 



^ i7rr72 2 j_ 

2^ 



X^(X,/i2 



Xg{X,ii' 



(19) 



in agreement with the leading-order contribution in Eq. (|^. 

The next step is to investigate quantitatively how rapidly the asymptotic result 
is attained in practice. We focus first on the matrix element part, Eq. (|T^. Figure 1 
shows the ratio of the left-hand side (the exact result) to the right-hand side (the 
asymptotic result), as a function of the rapidity gap A. Evidently, the two agree 
to better than 10% for A > 3.3 . This is not unexpected, since we are testing here 
the size of the 'power correction' terms of order e~^. We can also test the effective 
subprocess approximation, by comparing the qg — > qg and qq -^ qq amplitudes, 
scaled by 9/4 and (9/4)^ respectively, to the asymptotic gg — *• gg result. These ratios 
are shown as the dashed (qg) and dash-dotted (qq) lines in Fig. 1. The approach to 
the limiting form is very similar to the gg amplitude, indicating that in the large A 
region the effective subprocess approximation is valid. 

We have already discussed how when M and ^/s are fixed the bulk of the A 
dependence comes from the parton distributions. This is illustrated in Fig. 2, where 
the cross section of Eq. (0), scaled by M^, is shown as a function of A. We have 
chosen ^i = 1.8 TeV, and used the latest 'MRS(H)' partons with A^ = 230 MeV. 
The scales in the running coupling and in the parton distributions are both set equal 
to M. The solid curves correspond to the exact gg -^ gg cross section, for the 
representative values M = 10 GeV and M = 30 GeV. Also shown (dashed lines) 
are the asymptotic cross sections defined by Eq. (|1^). There is a broad range of A 
where the shape of the exact cross sections is reasonably well approximated. The 
normalization, however, is too high by a factor of order two.Q This can be traced to 
the fact that Eq. ([T8|) is only valid when both X is small and the function xg{x, /i) is 
slowly varying. If xg ~ a;"'*' in the relevant x region, then an error of order (1 + A) is 
made in the normalization. 

2.2 Comparison with exact O(a^) calculation for (j) ^ 

Expanding the exponential inside the integral in Eq. (^ allows us to read off the 
differential 0-distribution at next-to- lowest order, i.e. O(a^): 



(20) 

NLO 



_ (OLsCa^ 
V TT J 


)' YA/(0) 


1 +CO 

27r Jt-^ 


e^"^ c„ , 



, ^9 do" 
d0 

where 

/(</>) = :f E e^"'c„, (21) 



■^In the range 2 < A < 6 for Af = 10 GeV the ratio of the exact to the approximate cross section 
always hes in the range 0.45 to 0.55. 



and the Fourier coefficients are 

1 /■+°° dz 



-Xniz) = 2[V^(l)-^(l + i|n|)]. (22) 



" 27r 7-00 z^ + i 

The coefficients can be calculated explicitly by contour integration in the complex z 
plane: 

Cn — C-n , Cn+2 — C„ — — - [n > 0) , 

2 + n 
co = 0, ci = 4(log2-l). (23) 

The ffist point to note is that since cq = 0, the integral of /(</>) vanishes, i.e. 

M' a^Lo = . (24) 

To calculate the distribution we have to substitute the c„ coefficients into Eq. ( ^Ij) 
and sum the Fourier series. We have found it easier, however, to start from the 
solution to the BFKL equation in transverse momentum space, where for 7^ (see 
Appendix), 

M^ r°o /-oo dA;!,! dk^2 






2-71 JNP JNP kj.^kj.2{,kj.^ + k'?p2 — 2kTlkT2 COS0) 

= - [log(2(l - COS 0)) + (tt - 0) cot 0] for < < TT , (25) 

TT 

and f{—<f)) = /(0)- The singular behaviour / ~ (p^^ as — > 0, which corresponds 
to almost back-to-back jets, arises from soft emission of the third (gluon) jet ^. It 
is cancelled in the total cross section by a virtual gluon contribution proportional to 
S{(f)). This can be taken into account by invoking the standard 'plus prescription', 
i.e. /(0) ^ [/(0)]+ where 

d0 9{<l>) [fm+ = r # {9{<jy) - 9(0)) f{<f>) . (26) 

J —n 

It is important to note that the distribution we have derived is a leading loga- 
rithm result, in that it corresponds to retaining only the leading a^A contribution 
and ignoring corrections of order a^, a^e"^ etc. To study the validity of this approx- 
imation, we can compare the result in Eq. (|25|) with an exact calculation based on 
the complete 99 —>■ 999 matrix element. The analogue of the 2^2 cross section 
Eq. (IID is 



da 



dyidy2 



1 r ... r ... ri^ 



/ dk^i / dk^2 I I dys 



1. 5127r4 7m2 ..„- .-^. 

X [xi9{xi, fl"^) X2fi'(x2,/i^)] 



X s 



j:\M99 ^999)\' (27) 



where 






Xi 
S 



k'^i + k'^2 ~ 2A;ti^T2 cos0 
{kTie2^ + A;r2e~2^ + A;T3e^«)/v^ 

/c^i + A;|,2 + ^Ts + 2kj'ikT2 cosh A 

+2kTikT3 cosh(iA - t/s) + 2kT2kT3 cosh(iA + t/s) 



(2^ 



According to the BFKL analysis, this cross section should have the asymptotic limit 



da 



dyidy2 



yi= 



-y2-- 



4- 



a.CA 



IT 



TT 



2m2 



2M2 



A/(0)[X^(X,/i^)] 



(29) 



We can see how this behaviour arises: the matrix element in ( ^7\) is dominated by 
configurations where the third gluon jet is produced centrally, i.e. I1/3I <^ |A. With 
the matrix element approximated by its value at 2/3 = 0, the 2/3 integral gives the 
overall factor of A, and the remaining kxi integrals give the function /(0), as in 
Eq. (^). The parton distributions are again dominated by their values at Xi = X2 = 
X = 2Mcosh(|A)/A/s, as for the leading order cross section. 

We can study the approach to the asymptotic result in two stages. First, at the 
subprocess level, we can compare the cross section 



d(T 



dyidy2d(f) 



yi=-y2=2'^ 



51277^ 



Af2 



d4i dk^2 \ dy^s ^Y. \^(9g -^ 999)\' 



with the asymptotic form 

da 



dyidy2 



yi=~y2=-^^ 



dsCA^^ 



TT 



TT 



2M2 



A/(0) 



(30) 
(31) 



Figure 3 compares the distribution calculated from Eq. (|30D with the function 
/(0), for A = 4,8,12. The exact calculation has been scaled by the same factors 
multiplying /(0) on the right-hand side of Eq. (|3TD, so that the exact and approximate 
distributions should coincide in the limit A — > 00. The results confirm the approach 
to the asymptotic distribution. At small 0, the asymptotic behaviour is already a 
good approximation for A = 4, while the convergence is slower at large 0. This 
is presumably because at large the third gluon can have significant transverse 
momentum and energy, thus invalidating the approximations under which the BFKL 
equation is derived and leading to sizeable sub-asymptotic corrections. 
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Figure 4 makes the same comparison at the cross-section level, for pp collisions 
with ^/s = 1.8 TeV, M = 10 GeV and A = 4,6,8. The curves correspond to the 
exact and asymptotic cross sections of Eqs. (^) and (^) respectively. The constraints 
2^1; 2:2 < 1 now give upper limits on the transverse momentum integrals. As for the 
leading-order case, the normalization is overestimated by the asymptotic form, but 
evidently the shape of the distribution is reasonably well approximated even at 
moderate A. The rapid change in the distribution with increasing A supports our 
assertion that the azimuthal distribution of the jets should be a more reliable indicator 
of BFKL behaviour than the overall 0-integrated cross section. 

2.3 All-orders distribution 

Through next-to-lowest order, then, we have 
^2da(E,M) _ fasCA\^ n^ 



F(.^,A) = m + (^) A 1/(0)]+ + ... , (32) 

where the . . . represent terms 0((asA)") with n > 2. Formally, the first few terms 
in this series will be a good approximation to the all-orders distribution provided 
A ^ 1 and a^A -C 1. As the second of these inequalities is relaxed, higher order 
terms become more and more important. The inclusion of all terms of the form 
(a^A)", via Eqs. (|^J^), requires a numerical calculation and will be discussed below. 
First, we consider the limit a^A ^ 1, where an analytic approximation can again be 
obtained. 

To calculate the distribution in the asymptotic limit a^A — ^ cxd, we return to 
Eq. (^) and use a saddle-point method [Q] to evaluate the Fourier coefficients in the 
large t = agCA^ji^ limit. We expand the Xn(^) about the saddle point at z = 0, 



Xn{z) =an- hnZ^ + ... , (33) 



where 



2 

ao = 2 log 2, ai = 0, a„+2 = a„ - -- — , {n > 0) 

1 + n 

bo = 7C(3), 6i = C(3), 6„+2 = &n - ,^ ^ ,3 , {n > 0) (34) 

(1 + n)-^ 

from which the asymptotic t — > oo behaviour follows, 

Cr.it) ~ -j^ e^""* . (35) 

9 



Figure 5 shows the first seven C„ coefficients as functions of t. Note that because 
a„ < for n > 1, all but the Co coefficient tend to zero as t ^ oo. The asymptotic 
distribution is then obtained by substitution in Eq. (^), 



A) 



1 g41og2t ^ 2 COS. 



1 

2^ 



+ - g4(log2-l)t _j_ 



2 COS 
i7r"(7C(3) - 8)t 



(36) 



Asymptotically, then, the distribution becomes flat. The emission of an infinite 
number of soft gluons has completely smeared out the back-to-back correlation ex- 
hibited by the lowest contributions to the perturbation series. Note that we have also 
reproduced the asymptotic result for the 0-integrated cross section |l|, |^: 

M^a(„.,M^ A) = (^' ^ CM _ (^' ^ ^=4= e«'«- . (37) 

For large A, therefore, we obtain the result given in Eq. (|^) with 

(y 
A = ^4CAlog2 = 0.5, (38) 

TT 

for as = 0.19. The 0-integrated cross section was studied in some detail in Ref. 0, 
where an analytic approximation valid for t ^ 1 was derived. Figure 6 shows the 
function Co(t), (i) computed exactly using Eq. (|) (solid line), (ii) according to the 
analytic approximation of Ref. (dotted line), and (iii) in the asymptotic limit, 
Eq. (H i (dashed hue). 

We have so far obtained analytic approximations for the small t and large t be- 
haviours of the differential (p distribution. The distribution at arbitrary t requires 
a numerical calculation of the sum and integral in Eqs. (§,^. Thus, Fig. 7 shows 
the function F(0,A) defined in Eq. (D for t = a,C^A/7r = 0.25,0.5,1.0,1.5. We 
see very clearly the transition from a sharply peaked distribution at small t — recall 
that F = 6{(f)) at t = — to a larger, flatter distribution as t increases. As the 
rapidity gap widens, the emission of more and more soft gluons uniformly 'fills in' 
the distribution at large 0. 

3 Predictions for pp collisions at 1.8 TeV 

The most direct test of the BFKL perturbative pomeron behaviour is the rise in 
the subprocess cross section with increasing rapidity gap A, i.e. M'^a ~ exp(AA). 

10 



However, as discussed in the Introduction, one cannot yet regard this as a precision 
prediction of the theory. One particular issue concerns the inclusion of a running 
coupling in the BFKL analysis, i.e. as — > OsI^t)- This prevents the integrals being 
extended down to fc^ = 0, thereby inducing a weak dependence on an infra-red cut- 
off parameter, see for example Ref. 0. Increasing this cut-off reduces the phase 
space for the soft gluon emission and weakens the growth in the cross section with A. 
In addition, the subleading logarithmic corrections to the BFKL result are not yet 
known. As we have stressed, at fixed hadron collider energy the BFKL behaviour is 
anyway masked by additional dependence on A coming from the parton distributions. 
To investigate this latter effect quantitatively, we show in Fig. 8 the cross section of 
Eq. (|l|) at \fs = 1.8 TeV as a function of A for two choices of minimum jet transverse 
momentum, M = 10 GeV and M = 30 GeV. The dashed curves correspond to 
the leading order contribution to o", i.e. the first term on the right-hand side in 
Eq. (H), while the solid curves are the all-orders BFKL result, corresponding to a 
given in Eq. (^). The parton distributions are the latest MRS(H) set |P], with A^^ = 
230 MeV, which are consistent with both the recent HERA and the fixed-target F2 
measurements. Evidently, the shapes of the lowest-order and all-orders distributions 
are quite different. Notice that the x-dependence of the parton distributions more 
than compensates the BFKL rise in the subprocess cross section, so that the net effect 
is a cross section which decreases as a function of A. However, we should recall from 
Section ^]T| that the lowest-order approximation is only a good representation of the 
shape of the exact lowest order 2 — ;> 2 cross section for large A > 4. The signature for 
BFKL behaviour is therefore a slower fall-off of the cross section with increasing A 
than predicted by the leading order (exact or approximate) cross section. The size of 
the effect can be gauged from Fig. 8. Whether the difference is detectable in practice 
depends on the precision with which jets at large rapidity can be reconstructed and 
measured experimentally. 

We turn next to the distribution in the azimuthal angle difference of the two jets. 
Figures 9 and 10 show the = vr — (f)jj distributions, at fixed A, for ^/s = 1.8 TeV 
with (a) M = 10 GeV and (b) M = 30 GeV. In Fig. 9 the cross section itself is 
shown, while in Fig. 10 the distributions are normalized to have unit area for each 
value of A. The trend is that as A increases, the cross sections get smaller and the 
distributions become flatter in cj). The higher the transverse momentum cutoff M, the 
faster the decrease with A and the slower the approach to the flat distribution. For 
M = 10 GeV, values of A up to about 8 appear to be accessible, at least in principle. 
For this rapidity gap, the (f) distribution is almost fiat. A simpler representation of 
the flattening behaviour is provided by the average of cos 0, which is proportional 
to the Ci coefficient of Eq. (j^), i.e. 

(cos0) = ——r , t = (39) 

Co[t) TT 
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Figure 11 shows this average as a function of A, for M = 10 GeV and M = 30 GeV. 
The approach to flatness ((cos0) -^ 0) is slower for the higher cut-off because of the 
smaller coupling constant. 

4 Conclusions 

In this study we have focussed on two important aspects of BFKL perturbative 
pomeron behaviour. The first is the rise in the cross section as the rapidity gap 
between two moderate pt jets increases, as first discussed in Ref. 0. The difficulty 
with this signature is that in order to circumvent the additional dependence on ra- 
pidity induced by the parton distributions, it is necessary to scale up the collider 
energy as the rapidity gap increases. At fixed collider energy, the dominant effect at 
large rapidity gap is the suppression of the cross section by the fall-off in the parton 
distributions as x ^ 1. After allowing for the effects of jet reconstruction and mea- 
surement in the detector, it is not clear whether the relatively small effects of the 
BFKL behaviour can be observed. This is not, however, as severe a problem for the 
second important feature of the BFKL pomeron — the weakening of the correlation 
in the azimuthal angle of the two jets. The distribution in the azimuthal angle dif- 
ference changes from being back-to-back for central jet pairs with a small rapidity 
difference, to an asympototically fiat distribution as the jets separate in rapidity. 

We have presented predictions for the (j) distribution at different rapidity gaps A 
in pp collisions at \/s = 1.8 TeV, based on the fixed-coupling solutions of the BFKL 
equation. However, there are very likely non-negligible sub-leading corrections to this 
behaviour. We have investigated some of these by comparing the leading-logarithm 
predictions with those based on the exact low-order matrix elements. From this 
comparison, it seems that the leading behaviour may already be dominant for rapidity 
gaps as small as 4. 

Of course, we have not included such important effects as smearing of the jet ener- 
gies and angles by the jet algorithms used in the actual experiments. A more precise 
analysis would require a correspondong smearing of the distribution, which would 
inevitably weaken the correlations implied by perturbation theory alone. The size of 
this smearing could be estimated, for example, by comparing the actual distribution 
of two central jets — which should be well-described by lowest order perturbation 
theory — with the naive 5(0) expectation. If the smearing was parametrizable by, 
say, a gaussian distribution, this could be folded in to the Fourier coefficients in the 
perturbative predictions. We have not performed such an analysis here, since the 
form of the smearing is presumably detector and jet algorithm dependent. 

In summary, therefore, it would be interesting to measure the azimuthal angle 
distribution of the two-jet inclusive cross section as a function of the jet rapidity gap, 
to see if the data are at least qualitatively in line with the flattening of the distribution 
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predicted by the BFKL equation. In this study we have only considered the case oipp 
coUisions at 1.8 TeV. However, basically the same behaviour should also be manifest 
in any high-energy collider with quarks and gluons in the initial state. In particular, 
photoproduction of jet pairs at HERA could also be a useful place to look for evidence 
of the BFKL behaviour, and, of course, high-energy proton-proton colliders such as 
the LHC will allow a much larger range of rapidities to be covered. 



Appendix: Solution of the BFKL equation 

In this Appendix, we present a brief derivation of the solution of the BFKL equation 
for the two-jet inclusive cross section. More details can be found in Refs. fl], H). 

We start from the subprocess cross section gg — » gg + X, where the two final 
state gluons are produced with transverse momenta kTi at large rapidity separation 
A, and X represents additional soft gluons. The differential cross section can be 
written, following the notation of |H], 



do- _ ^2^.2 /(fc Tl,fcT2, A) 






The Laplace transform of the function /, 

f{kTi,kT2,uj)= dAe--^f{kTukT2,A), (A2) 

Jo 

satisfies the BFKL equation [Q] 



ujf{kTi,kT2,uj) = 6{k^i-k'^2)^{<Pi-<p2)+ {^^ 



vr^ 



X 



d'k. 



{kri - kr) 



2 



f{kT,kT2,Uj) 



kTlfikTl,kT2,^) 

kl + {kri - krf 



(A3) 



This equation can be solved in closed form by introducing the Fourier transform of 
/ with respect to 0i — 02 and \og{k'^i/k'^^: 

fikru kT2, ^) = 7^ E e*"^'^^-^^) 7^ r d-^ e-^^i"^^'^"/^-) /„(^, u). (A4) 

Substituting this into Eq. ( |A3| ) gives 

ujfn{z,uj) = (^kTik'^2) ^ + u}o{n,z) fn{z,uj) , (A5) 
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where 

u,{n,z) ={^) 2xn{z), (A6) 



with the function Xn given in Eq. (^. Performing the inverse transform of Eq. (|A2| ) 
then gives 

1 1 /"oo 1 1 

/(A;n, fcT2, A) = - E e^"^ TT / ^^ (4i)~2-^^(42)"2+*^ e2*'^"(^) , (A7) 

ZTX „ ZTT J-oo 

where = 0i — 02 and t = asC^A/vr. The final step is to integrate the transverse 
momenta over the range M^ < k'^^ < oo, 

jr roo /"OO rn Apr 

2 Jnp Jap j~ 



2Jnp ''Jap '"7-^ ^ d'^kTid^kT2 
alCl-K r , . I ^^ .-...A. 1 r^ dz 



-■A' 



I>i^^-"'rJl^^"'^'"' (*^ 



2M2 7-^ 27r ^ 27r J-oc z'^ + ^ 



which gives the subprocess cross section of Eqs. (§,|,|]). 
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Figure Captions 

[1] Ratio of the exact 2 — > 2 subprocess cross sections defined in the text (Eq. ([17|) ) 
to the asymptotic gg -^ gg scaling form given in Eq. (^, as a function of the 
rapidity gap A. The curves are (i) gg — > gg (solid line), (ii) qg — > qg multiplied 
by 9/4 (dashed line), and (iii) qq — > qq multiplied by (9/4)^ (dot-dashed line). 

[2] The lowest order 2 — i> 2 cross section of Eq. (|I3]) as a function of the rapidity 
gap A, at yi = 1.8 TeV. The MRS(H) parton distributions [§ with A^ = 
230 MeV are used. The solid curves correspond to the exact gg -^ gg cross 
section with minimum jet transverse momenta M = 10 GeV and M = 30 GeV. 
Also shown (dashed lines) are the asymptotic cross sections defined by Eq. (p!9|). 

[3] The asymptotic azimuthal angle distribution /(</>) of Eq. (pSf), compared to the 
distribution calculated from the exact gg -^ ggg matrix element for A = 4, 8, 12 
(dashed lines). 

[4] The dependence of the differential cross section at 0{al) on the azimuthal 
angle difference of the two jets with rapidity gap A = 4, 6, 8, in pp collisions 
at yi = 1.8 TeV with M = 10 GeV. The MRS(H) parton distributions g 
with Ajj^ = 230 MeV are used. The solid curves correspond to the exact cross 
section (Eq. (0)), and the dashed curves to the asymptotic approximation of 
Eq. m. 



[5] The Fourier coefficients C„ (n < 6) defined in Eq. (||) as functions of t = 

asC^A/TT. 

[6] The Fourier coefficient Co(t), which gives the 0- integrated cross section, (i) 
computed exactly using Eq. (||) (solid line), (ii) according to the analytic ap- 
proximation of Ref. P (dotted line), and (in) in the asymptotic limit, Eq. ( pT)) 
(dashed line). 

[7] The azimuthal angular distribution function F(0, A), defined in Eq. (^, for 
t = asCAA/n = 0.25, 0.5, 1.0, 1.5. 
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[8] The cross section dcr/d|/id|/2(|/i = —y2 = |A) of Eq. (|I|) at ^/s = 1.8 TeV 
as a function of A, for two choices of minimum jet transverse momentum, 
M = 10 GeV and M = 30 GeV. The dashed curves correspond to the leading 
order contribution to a, i.e. the first term on the right-hand side in Eq. (0), 
and the sohd curves are the all-orders BFKL result, corresponding to a given 
in Eq. (jS^). The MRS(H) parton distributions § with a|^ = 230 MeV are 
used. 

[9] The = TT — (f)jj distributions, at fixed A, for y/s = 1.8 TeV with (a) M = 
10 GeV and (b) M = 30 GeV. The MRS(H) parton distributions |[ with 
A^ = 230 MeV are used. 

[10] As for Fig. 9, but with the distributions normalized to unit area at each A. 

[11] The average azimuthal angle difference (cose/)) = — {cos (pjj) as a function of A, 
for M = 10 GeV and M = 30 GeV. 
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